arXivil506.02883vl [astro-ph.GA] 9Jun2015 


Astronomy & Astrophysics manuscript no. document ©ESO 2015 

June 10, 2015 


Infrared-faint radio sources remain undetected at far-infrared 

waveiengths^ 

Deep photometric observations using the Herschel Space Observatory 

A. Herzog^’^’^, R. R Norris^, E. Middelberg\ L. R. Spitler^’^, C. Leipski^, and Q. A. Parker^’ 

^ Astronomisches Institut, Ruhr-Universitat Bochum, Universitatsstr. 150, 44801 Bochum, Germany 
e-mail: herzog@astro. rub. de 
^ Macquarie University, Sydney, NSW 2109, Australia 

^ CSIRO Astronomy and Space Science, Marsfield, PO Box 76, Epping, NSW 1710, Australia 
^ Australian Astronomical Observatory, PO Box 915, North Ryde, NSW 1670, Australia 
^ Max-Planck-Institut fur Astronomie, Konigsstuhl 17, 69117 Heidelberg, Germany 

^ Department of Physics, Chong Yeut Ming Physics Building, The University of Hong Kong, Pokfulam, Hong Kong 
Received 26 November 2014 / Accepted 27 May 2015 


ABSTRACT 


Context. Showing 1.4 GHz flux densities in the range of a few to a few tens of mJy, infrared-faint radio sources (IFRS) are a type of 
galaxy characterised by faint or absent near-infrared counterparts and consequently extreme radio-to-infrared flux density ratios up to 
several thousand. Recent studies showed that IFRS are radio-loud active galactic nuclei (AGNs) at redshifts > 2, potentially linked to 
high-redshift radio galaxies (HzRGs). 

Aims. This work explores the far-infrared emission of IFRS, providing crucial information on the star forming and AGN activity of 
IFRS. Furthermore, the data will enable examining the putative relationship between IFRS and HzRGs and testing whether IFRS are 
more distant or fainter siblings of these massive galaxies. 

Methods. A sample of six IFRS was observed with the Herschel Space Observatory between 100 pm and 500 pm. Using these results, 
we constrained the nature of IFRS by modelling their broad-band spectral energy distribution (SED). Furthermore, we set an upper 
limit on their infrared SED and decomposed their emission into contributions from an AGN and from star forming activity. 

Results. All six observed IFRS were undetected in all flve Herschel far-infrared channels (stacking limits: cr = 0.74 mJy at 100pm, 
cr = 3.45 mJy at 500pm). Based on our SED modelling, we ruled out the following objects to explain the photometric characteristics 
of IFRS: (a) known radio-loud quasars and compact steep-spectrum sources at any redshift, (b) starburst galaxies with and without an 
AGN and Seyfert galaxies at any redshift, even if the templates were modifled, and (c) known HzRGs at z < 10.5. We And that the 
IFRS analysed in this work can only be explained by objects that fulAl the selection criteria of HzRGs. More precisely, IFRS could be 
(a) known HzRGs at very high redshifts (z > 10.5), (b) low-luminosity siblings of HzRGs with additional dust obscuration at lower 
redshifts, (c) scaled or unsealed versions of Cygnus A at any redshift, and (d) scaled and dust-obscured radio-loud quasars or compact 
steep spectrum sources. We estimated upper limits on the infrared luminosity, the black hole accretion rate, and the star formation rate 
of IFRS, which all agreed with corresponding numbers of HzRGs. 

Key words. Techniques: photometric - Galaxies: active - Galaxies: high-redshift - Galaxies: star formation - Infrared: galaxies 


1. Introduction 


A new class of extreme radio galaxies, which are characterised 
by their infrared (IR) faintness and their enormous radio-to-IR 
flux density ratios, has recently been found in deep surveys: 
the cl ass of infrared-faint radio sources (IFRS; [Norris et al.| 
,2006| ). While various suggestions had been presented to explain 
the extreme characteristics of these objects, recent observations 
clearly indicate that IFRS are radio-loud (RL) active galactic nu¬ 
clei (AGNs) at high redshifts. 

Infrared-faint radio sources were first discovered in the 
Chandra Deep Field South (CDFS) and the European Large Area 
Infrared space observatory Survey South 1 (ELAIS-Sl) fi eld of 
the Australia Telescope Large Area Survey (ATLAS) by Nor- 


* Herschel is an ESA space o bservatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 


ris et al.j ( |2006| ) and [Middelberg et al.j p008a| ), respectively. 


The radio maps at 1.4 GHz provided detections with fiux den¬ 
sities in the order of tenths to tens of mJy, whilst the co-l ocated 
Spitzer Wide-area Infrared Extragalactic Survey (SWIRE; [Lons 


[dale et aL]|2003| ) showed only a faint or, in most c ases, no IR 
counte rpart at 5.6pm with a noise a = 1 pJy. Later, jZmn et al. 
(2011|) defined the class of IFRS by two selection criteria: 


(i) radio-to-IR fiux density ratio S iaguzIS 3 , 6 jdm > 500, and 

(ii) 3.6yum fiux density 5 3 . 6 yum < 30yuJy . 


The enormous radio-to-IR fiux density ratios imply that 
IFRS are clear outliers from the radio-IR correlation. The second 
criterion implicates a distance selection, preventing ordinary ob¬ 
jects of redshift <1.4 from being included in the class of IFRS. 
In addition to the detections by [Norris et al.| ( |2006| ) and [Midd^ 


[berg et aL] ( [2008a[ ), IFRS were later found in the Spitzer extra- 


galactic First Look Survey (xFLS) field by Garn & Alexander 
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([2008|), in the Co smological Evolution Survey (COSMOS) by 
Zinn et aL] ( |2Qll| ), in the European Large Area IR space obser¬ 


vatory Survey North 1 (ELAIS-N p by Banfield et aLjpQllI ), 
and in the Lockman Hole field by |Maini et al.| ( |2013| ). In total, 
around 100 lERS have been found in these deep fields at a sky 
density of a few per square degree. 

Recently, |Collier et al. p014| ) found 1317 lERS in shal¬ 
low all-sky catalogues, that satisfied both selection criteria given 


above. They used the Unified Radio Catalog (URC; Kimball & 
Ivezic|2008|) based on the NRAO VLA Sky Suvey (NVSS;|Con 


don et al.||1998| ) and dat a from the all-sky W i de-Eield Infra red 


lERS are hotspots of spatially separated radio galaxies. If the 
IR detections were real counterparts of the radio detections, the 
redshift distribution of these objects would be expected to fol¬ 
low that of Eanaroff Riley type 2 (ERII) galaxies which peaks 
at z ~ 1. In contrast, if the IR detections were spurious asso¬ 
ciations to the radio emission, their redshift distribution would 
be expected to follow that of star forming galaxies, peaking at 


Survey Explorer (WISE; Wright et al.pOlO) . |Collier et al.| re- 
placed 53 . 6/^111 by S' 3 . 4 All lERS in that sample provide an IR 
counterpart at 3.4 pm and have radio flux densities above 7.5 mJy 
at 1.4 GHz, some of them exceeding 100 mJy. By this, the lERS 
in this sample are on average radio-brighter than the lERS found 
in the deep fields by around one order of magn itude. Since the 
WISE sensitivity is non-uniform, [Collier et al.| could only set a 
lower limit on the sky density of lERS with a radio fiux density 
above 7.5 mJy of ~ 0.1 deg“^. 

Eirst attempts to explain the existence of lERS included 
lobes or hotspots of radio galaxies, obscured star forming galax¬ 
ies (SEGs), pulsars or high-redshift RE galaxies. Later results 
(for an overview, see [Collier et al.|2014] and references therein) 
provided evidence that lERS are indeed RE AGN at high red- 
shifts. 

[Norris et ar]p007[ ) and [Middelberg et aL] ( [2008b| ) used very 
long baseline interferometry (VLBI) to confirm the presence of 
active nuclei in two lERS. Recently, [Herzog, A. et al.[ ( [2015[ ) de¬ 
tected 35 out of 57 lERS in VLBI observations, unquestionably 
confirming that at least a significant fraction of all lERS contains 
AGN. 

[Garn & Alexander| ( [2008| ) and [Huynh et~ar] ( [2010| ) present 
the first spectral energy distribution (SED) modelling of lERS 
and find that 3C sources (for example 3C 273) can reprodu ce the 
observational data when redshifted to z > 2. Eurthermore, [Garn[ 
[& Alexander[ and [Huynh et aL] conclude that obscured SEGs can¬ 
not explain the characteristics of lERS since the radio-to-IR fiux 
density ratios of lERS significantly exceed those of SEGs. This 
finding suggests the presence of an AGN in lERS, too. 

Recently, the first spectroscopic redshifts of lERS were 
found, confirmi ng the suggested high -redshift character of this 
class of objects. [Collier et ar] ( [2014| ) find 19 redshifts for their 
all-sky sample of lERS, all but one in the range 2 < z < 3. 
The outlier is expecte d to be a misidentifica tion or an SEG with 
an AGN in its centre. [Herzog et al.[ ( [2014[ ) measured redshifts 
of three lERS in the deep ATLAS fields based on optical spec¬ 
troscopy from the Very Large Telescope (VET) and found them 
to be in the range 1.8 < z < 2.8. The lERS with spectroscopic 
redshifts, both from [Collier et al.[ and [Herzog et al.[ lie at the IR- 
bright boundary of this class of objects with IR flux densities be- 


z ~ 0. Both cases are in clear disagreement with the observed 
redshift distribution of lERS, providing redshifts only in the 
regime 1.8 < z < 3.0. 

Based on the redshifts measured in their work, [Herzog 
[et al.[ p014[ ) present the first redshift-based SED modelling of 
lERS. They find the templates of an RL AGN (3C273) and 
a compact-steep spectrum (CSS) source (3C48) in agreement 
with all available photometric data, while an ultraluminous IR 
galaxy (ULIRG, Arp 220) or a Seyfert galaxy (Mrk231) are 
clearly ruled out. 

[Huynh et ar| ( [2010[ ) and [Middelberg et aL[ ( [2011| ) first pointed 
out a potential similarity between lERS and high-redshift radio 
galaxies (HzRGs). HzRGs are a class of powerful radio galax¬ 
ies (L 3 GHZ > 10^^ WHz“^) at high redshifts (z > 1). They are 
amongst the most massive and most luminous galaxies in the 
early universe and are expected to be the progenitors of the most 
massive galaxies in the local universe (e.g. [Seymour et al.|2007t 
[De Breuck et al.[[2010[ ). Moreover, HzRGs are known to host 
powerful AGNs and to show high star formation rates (SERs). 
Therefore, they are an important tool to study galaxy evolution 
and the interplay between AGNs and star forming activity. How¬ 
ever, HzRGs are rare objects; only around 100 of these objects 
are known in the entire sky, located at redshifts 1 < z < 5. 


Huynh et al. note that lERS have similar mid-IR to radio flux 


tween 14 p]y and 30 pJy at 3.6 pm or 3.4 pm. Both Collier et al 
and [Herzog et al.[ suggest that the lERS with measured spectro¬ 
scopic redshifts are the closest ones and that IR-fainter lERS are 
at even higher redshifts. 

Additionally, the results by [Collier et al.[ ( [2QT^ and [Herzog 
et al.[ ( [2014| ) finally disposed of the alternative hypothesis that 


density ratios to those of HzRGs. [Middelberg et al.[ studied the 
radio properties of lERS in the ELAIS-Sl field and find steep 
radio spectra between 2.3 GHz and 8.4 GHz with a median of 
a = -l.^and no index larger than -0.7. This is steeper than 
the general radio souce population (a = -0.86) and the AGN 
source population (a = -0.82) in that field. [Middelberg et al. 
find that the radi o spectra are even steep er than those of the sam¬ 
ple of HzRGs by [Seymour et al.[ ( [2007[ ), showing a median radio 
spectral index a = -1.02. HzRGs are the only objects known at 
significant redshifts sharing the extreme radio-to-IR flux den¬ 
sities of lERS in the order of thousands as shown by [Norris 
[et al.[ ( [2QTT] ). Eurthermore, [Norris et al.[ suggest that lERS might 
follow the correlation between 3. 6pm flux densit y and redshift 
found for the sample of HzRGs by [Seymour et al.[ similar to the 
K - z relation for other radio galaxies ( [Willott et al.|2003[ ). 

[Collier et~ar] ( [2014| ) and [Herzog et al. ( 2014 ) tested this hy¬ 
pothesis in the reachable redshift range between 2 and 3 and find 
lERS follow this correlation. This adds evidence to the sugges¬ 
tion that IR-fainter lERS are located at even higher redshift, po¬ 
tentially reaching z ~ 5 or 6. Eurthermore, [Herzog et al.[ show 
that the radio luminosities of lERS are in the same range as for 
HzRGs, although lERS lie at the lower bound of the radio lumi¬ 
nosity distribution of HzRGs. 

Besides these similarities between lERS and HzRGs, there 
is one significant difference. While HzRGs are rare objects of 
which only around 100 are known in the entire sky, lERS are 
much more abundant. The observ ed sky density of l ERS in the 
ATLAS fields is around 5 deg“ 


Zinn et al. 


( |2011| estimate a 

survey-independent sky density of (30.8 + 15.0) deg“^. 

Summarising, there is growing evidence that lERS are re¬ 
lated to HzRGs in the sense that lERS are fainter, but much 
more abundant siblings of HzRGs, potentially at even higher red¬ 
shifts. HzRGs are known to be vigorously forming stars while 
harbouring AGNs in their centres. Both components dominate 
the emission of galaxies in the IR regime. Observations at far- 
infrared (EIR) wavelengths are a key test for the hypothesis that 
lERS are related to HzRGs. 


The spectral index is defined as 5 oc y® 
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A. Herzog et aL: Infrared-faint radio sources remain undetected at far-infrared wavelengths 


Here, we analyse the IR and particularly the FIR regime of 
IFRS based on Herschel observations. Results in this wavelength 
range will enable us to further study the suggested link between 
IFRS and HzRGs. 

This paper is organised as follows. In Sect.j^ we describe the 
sample selection of the observed IFRS and the observations with 
Herschel. The calibration of the resulting data and the mapping 
is presented in Sect. In Sect. we describe the photometry 
and the consequent stacking analysis. We use the resulting flux 
density upper limits in Sect. for two different approaches to 
model. First, we model the broad-band SED of IFRS based on 


SED templates of known galaxies (Sect. |5.I| ). In a second ap¬ 
proach, we set an upper limit on the IR SED of IFRS in order to 
constrain their IR emission (Sect. |5.2| ). Based on the latter fitting, 
we estimate the black hole accretion rate and upper limits for the 
star formation rate. Finally, in Sect. we present our conclu¬ 
sion. In this paper, we use flat ACDM cosmological parameters 
Qa = 0-7, Dm = 0.3, Hq = 70 km s“^ Mpc“^ and the calculator 
by |Wright| ( |2^ . 


2. Sample and observations 

Six IFRS from the sample compiled by |Zinn et ak] ( |20II[ ) were 
photometrically observed with the ESA Herschel Space Obser¬ 
vatory ( [Pilbratt et al.||20I0| ), using the instrumen ts PACS (Pho- 
todetecting Array Camera and Spectrometer; |Griflin et al 


20I0|) and SPIRE (Sp ectral and Photometric Imaging Receiver; 


Poglitsch et al.|2010 


The sources were selected to be compar¬ 
atively bright in the radio regime, showing 1.4 GHz flux den¬ 
sities between 7 mJy and 26 mJy, and to provide high radio-to- 
IR flux density ratios (see Table also in the following). Since 
all observed sources are undetected in the SWIRE data between 
3.6//m and 24yum, their 5'i.4GHz/6'3.6yum lower limits based on 
SWIRE range between 2300 and 8700. Five observed IFRS are 
located in the field of ELAIS-SI and one in the CDFS. No IFRS 
located in the ELAIS-SI field has been detected in the X-ray 
XMM-Newton survey (flux limit of ~ 5.5 x 10“^^ erg cm“^ s“^ in 


the 0.5-2 keV band; [Puccetti et al.|2006| ) as mentioned by |Zinn| 
|et al.| ( |20lT] ). IFRS S703, located in the CDFS fleld, has not been 
covered by the CDFS Chandra 4 Ms survey ( |Xue et aL]|20II|) 


No optical counterpart is known for the six observed IFRS (Nor- 
ris et al.|2006[|Middelberg et al.|2008a| ). By these characteristics. 


the six obser ved sources are pro totypical for the class of IFRS. 
We note that |Maini et aI.| ( [20T^ And a 3 cr counterpart at 3.6 jum 
based in the Spitzer Extragalactic Representative Volume Sur¬ 
vey (SERVS; pal^duit et al.||20T^ for lERS S427 after the ob¬ 
servations with Herschel were carried out. 

Eigure shows the SWIRE 3.6yum maps of all lERS ob¬ 
served with Herschel, overlayed by the 1.4 GHz radio contours. 
lERS S427, S509, and S798 are point-like sources and show al¬ 
most no extended structure. Eor S703, in contrast, it might be 
suggested because of the additional radio component at around 
RA 03h35m33.5s Dec -27°26'40” and the bridge to the lERS it¬ 
self that the lERS is part of a double lobe radio source with the 
lERS and the additional component as radio lobes and the host 
galaxy in between those two components. However, the addi¬ 
tional radio component shows a weak 3.6pm counterpart, ex¬ 
cluding this component as a radio lobe. Therefore, S703 can be 
considered as a proper lERS. 

Infrared-faint radio source S973 
tended structure and has been listed 
nents (RA00h38m44.723m Dec-43° 

RA 00h38m43.489s Dec -43° I9'26.94") 


an 


ex¬ 


shows 

as two compo- 
I9H4. 58^^ and 
by [Middelber^ 
et al.| ([2008a|), with the source position in between the two 




0.05 0.1 0.15 0.2 

3.6 flux density [MJy sr ^ ] 


Fig. 1. SWIRE 3.6pm maps (greyscale; Lonsdale et al. 

2003 

1 of all 

six IFRS observed with Herschel overlaved bv the 1.4 Gf 

z rad 

0 con- 

tours (ATLAS data release 3; cr ^ 20yuJy beam 

Franzen et al. 

2015I. 


Shown are the contours at 2cr, 4cr, 8cr, 16cr, etc. Top down and from 
left to right: S703, S427, S509, S749, S798, and S973. 


components. Here, however, we considered the position of the 
second, radio-brighter component to be the centre of the host 
with the weaker component more likely to be a jet. However, the 
choice of host position does not affect our overall conclusion. 
The position of the brighter component is listed in Table 

Infrared-faint radio source S749 shows an extended struc¬ 
ture in the radio map, too. However, this source is located close 
to the edge of the ATLAS ELAIS-SI fleld and, therefore, the 
noise is higher than for the other lERS discussed in this work. 
In the scheduling of the Herschel observations, a slightly wrong 
position of lERS S749 had been used. Thus, this source is not 
located in the centre of the Herschel maps but closer to the edge, 
resulting in a higher noise compared to the other sources, par¬ 
ticularly at shorter wavelengths. Therefore, we will exclude this 
source from our stacking analysis in Sect. 

The PACS observations of these six objects were carried out 
in July and December 2011 (observation IDs 1342224373, 
1342224374, I3422336I5, I3422336I6, I3422336I7, 

I3422336I8, I3422336I9, 1342233620, I34223362I, 

1342233622, 1342233623, 1342233624). Using the mode 
mini-scan map, each observation was divided into two parts, 
arising from two different scan angles. Both scanning parts were 
centred on the source, observing simultaneously at lOOyum and 
160 pm with a total on-source time of 22.5 minutes for each 
source. 

The SPIRE observations were carried out in December 2011 
and January 2012 (observation IDs 1342234729, 1342234730, 
I34223473I, 1342234732, 1342234733, 1342238290). Using 


Article number, page 3 of 


15 













































































































































ASlA proofs: manuscript no. document 


Table 1. Sample of IFRS observed with Herschel. 


lERS 

ID 

RA 

J2000.0 

<^RA 

[arcsec] 

Dec 

J2000.0 

t^Dec 

[arcsec] 

s 1,4 GHz 

[mJy] 

^ 3.6^m 

IpJy] 

‘^24jum 

iMJy] 

>5 1.4GHz/5'3.6/im 

references 

S703 

03:35:31.025 

0.11 

-27:27:02.20 

0.12 

26.1 

< 2.04 

< 115 

> 12784 

(1), (4) 

S427 

00:34:11.592 

0.01 

-43:58:17.04 

0.00 

21.4 

1.77 + 0.54 

< 150 

12000 

(2), (4) 

S509 

00:31:38.633 

0.01 

-43:52:20.80 

0.01 

22.2 

< 3 

< 150 

> 7400 

(2), (3) 

S749 

00:29:05.229 

0.04 

-43:34:03.94 

0.04 

7.0 

< 3 

< 150 

> 2337 

(2), (3) 

S798 

00:39:07.934 

0.04 

-43:32:05.83 

0.03 

7.8 

< 3 

< 150 

> 2597 

(2), (3) 

S973 

00:38:43.489 

0.04 

-43:19:26.94 

0.04 

9.1 

< 2.70 

< 150 

> 3385 

(2), (4) 


Notes. Infrared-faint radio source S703 is located in CDFS while the other five sources are located in the E LAIS-Sl field. All six IFRS are 
undetected at 3.6//m in the SWIRE survey (3cr ~ 3 piy) as presented by |Zinn et al.| { |20Il) . |Maini et al.| ( |20I3| ) analysed the deeper SERVS data 
for three lERS observed with Herschel and find one week counterpart. All upper limits represe nt 3cr. S973 consis ts of two radio components. 
Here, we used the position of the strongest component as the position of the source, in contrast to |Middelberg et al!] who used the centre between 
both components. Positions, position uncertainties, radio and 24 pm fiux densities were taken from the first reference listed in each row, S 3 . 6 //m and 
S i. 4 GHz/‘^ 3 . 6 ;im from the second reference. 

References. (l) |Norris et ari ( [2006t ; (2) |Middelberg et al.| ( [2008al ); (3) |Zinn et al.| ( [20TT] ); (4) |Maini et 


the observing mode SpirePhotoSmallScan, the sources were 
observed simultaneously at 250yum, 350yum and 500pm for 
12 minutes each. 

We recognise that, since these observations were made and 
this analysis was performed, additional Herschel data on IFRS 
have been released into the public domain by the Herschel Multi- 
Tiered Extragalactic Survey (He rMES; [Oliver et al.||2012| ), the 
PACS Evolutionary Probe (PEP; [Lutz et al.|2Qll| ), and the Cos¬ 
mic Assembly Near-infrared Deep Extragalactic Legacy Sur¬ 
vey (CANDLES ; [Grogin et al.|20TT] ). These data have not been 
used for the study presented in this work, but will be discussed 
in a future paper. 


3. Data calibration and mapping 

The data sets were calibrated using the Herschel Interactive Pro¬ 
cessing Environment (HIPE; [Ott|2Q10 version 12.1.0). We fol¬ 
lowed the steps presented in the PACS photometer pipeline for 
deep survey maps provided within HIPE to process the PACS 
data sets from level 0 to level 2. During the processing to level 1, 
a mask was created and applied in the task highpass filter 
to prevent artefacts in this median-subtraction process arising 
from nearby bright sources. The position of the lERS has also 
been masked. We obtained the best results using a highpass 
filter radius of 15". The final maps were built using the task 
photProject, setting the parameter pixfrac to 0.5, both for 
100yum and 160/im. Output pixel sizes of 1.5" (for 100yum) and 
2 .1" (for 160yum), respectively, were found to provide the best 
maps. Both quantities are parameters of the task photProj ect 
which is based on the Drizzle method ( [Eruchter & Hook|2002[ ). 
Einally, we matched the two maps of each object, obtained from 
the two different scan angles, using the task mosaic. The final 
PACS maps have a size of around 3.5' x 6.5'. The central re¬ 
gion with a diameter of around 50"provides the highest and an 
almost uniform sensitivity. Eigure (upper subplots) shows the 
final PACS maps of lERS S509. 

Correspondingly, for the SPIRE data sets, we followed the 
steps in the appropriate HIPE standard pipeline—the “Photome¬ 
ter Small Map user pipeline”—to calibrate the data. The maps 
were built using the task Destriper, an iterative algorithm to 
remove the baseline signal in the timeline data. In this task, the 
implemented deglitcher was used. We used the algorithm naive 
mapping implemented in HIPE to build the final maps, project¬ 
ing the full power seen by each bolometer timeline step to the 
closest sky map pixel with the pixel sizes given in Table The 
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Fig. 2. Final Herschel maps (greyscale) of IFRS S509 overlayed by the 
1.4 GHz radio conto urs (ATLAS data release 3; cr ^ 20yL/Jy beam"^; 
[Franzen et al.|2015] l. Shown are the contours at 2cr, 4cr, 8cr, 16cr, etc. 
Top down and from left to right: 100 pm, 160 pm (both PACS), 250 pm, 
350yum, and 500yum (SPIRE) maps. None of the five maps provides a 
detection of S509. 


usable area of the final SPIRE maps is around 8' x 10'. Eig¬ 
ure (lower three subplots) shows the resulting SPIRE maps of 
lERS S509. 

























































































































A. Herzog et aL: Infrared-faint radio sources remain undetected at far-infrared wavelengths 


Table 2. Characteristics of the instruments and parameters of the data reduction and photometry, used in HIPE for the different channels. 


Band 

FWHM 

[arcsec^] 

Pixel size 
[arcsec pixel" 

Exposure time 
[s] 

PACS lOOyum 

6.69 X 6.89 

1.5 

1354 

PACS 160yum 

10.65 X 12.13 

2.1 

1354 

SPIRE 250yum 

17.6 X 17.6 

6 

721 

SPIRE 350yum 

23.9 X 23.9 

10 

721 

SPIRE 500 yum 

35.2 X 35.2 

14 

721 


Notes. The numbers of the FWHM are taken from ‘TAGS Observer’s Manual” version 2.5.1 and from “The SPIRE Handbook” version 2.5, 
respectively. 


4. Photometry and stacking 

4.1. Photometry 

We used the HIPE task AnnularSkyAperturePhotometry for 
photometry on the PACS maps, performing aperture photometry 
on a chosen target based on a circular aperture. The related task 
PhotApertureCorrectionPointSource corrected the mea¬ 
sured flux density for the finite size of the applied aperture, yield¬ 
ing the requested flux density of the nominal source. We used 
aperture radii of 5.6"and 10.5"for the PACS maps at 100/im 
and 160/im, respectively. The background was estimated und 
subtracted based on a ring between 20" and 25", and 24" and 
28", respectively. The flux density uncertainties were obtained 
from 100 similar apertures randomly placed on the map where 
the coverage was at least 75% of the maximum coverage. We 
fit a Gaussian to the histogram of these background flux densi¬ 
ties, providing the Gaussian width cr. We repeated this procedure 
ten times and took the median of these Gaussian widths as the 
uncertainty of the measured flux density. For S749, where the 
IFRS was not in the centre of the Herschel maps as discussed in 
Sect. 1^ we placed 100 similar apertures in a field with a max¬ 
imum distance to the source of five times the full width at half 
maximum (FWHM) of a point source to estimate the flux uncer¬ 
tainty. 

Our six nominal objects were not detected at lOOyum or at 
160yum. Table [^summarises the resulting point source flux den¬ 
sity uncertainties cr for the six observed IFRS. 

Photometry on the SPIRE maps was carried out using the 
HIPE task sourceExtractorSussextractor based on the 
SUSSEXtractor algorithm ( [Savage & Oliver | [2007 1 ). Adapting 
this task to the maps, sources were extracted, providing data 
about brightnesses, positions, and related errors of the sources 
as well as a cleaned map. The extracted sources were checked 
for their angular distance to the known position of IFRS (< 
0.5 FWHM). We repeated this procedure as a second level source 
extraction on the cleaned map. Again, we checked the extracted 
sources for their distance to the IFRS’ positions. By this, we 
were able to extract sources which had been overseen in the first 
extraction step because of bright and close by sources in the map. 
These sources had been eliminated in the cleaned map. 

For five out of the six IFRS, the source extraction did not 
provide any FIR component in any SPIRE channel, either by the 
source extraction performed on the map itself or on the cleaned 
map. For S973, we found Scr detections in the cleaned maps at 
350 yum and 500 yum. The positions of these two FIR counterparts 
were in agreement with each other and were 12" distant from the 
IFRS radio position, which is only slightly below our distance 
criterion. However, we found that these FIR detections were very 
close to one source which was detected with a flux density of 
6.1 yuJy at 3.6yum in the SERVS survey. Because of those three 


reasons—weak flux, relatively large distance to the IFRS posi¬ 
tion, and overlap with a SERVS source—we suggest that these 
FIR counterparts are associated with the SERVS source and not 
with our SERVS-undetected IFRS. 

We summarise that none of the six observed IFRS was de¬ 
tected in any SPIRE channel. The flux density uncertainties were 
obtained from the cleaned map by fitting a Gaussian to the pixel 
values within a square of the size of eight times the FWHM of 
a point source (see numbers in Table [^. The SPIRE flux density 
uncertainties at the positions of the IFRS are summarised in Ta¬ 
ble O Our mean noise is lower than the overall confusion noise 
of 5.8mJy, 6.3 mJy, and6.8mJy at250yum, 350yum, and 500yum, 
respectively ( [Nguyen et al.||20I0| ). However, [Nguyen et al.[ also 
showed that the residual confusion noise after removing bright 
sources is lower. Our noises are in agreement with these num¬ 
bers. 


4.2. Stacking of Herschel maps 

Since no counterpart of an IFRS was detected in the Herschel 
maps, we performed a stacking analysis to search for a poten¬ 
tially weak counterpart slightly below the detection limit. The 
positional uncertainties for our sources are at least one order of 
magnitude lower than the pixel size of the Herschel maps (see 
Tables [^ an d|g. For each observed wavelength, we stacked the 
maps centred on the known position of the IFRS. However, we 
excluded IFRS S749 from the stacking analysis because of the 
higher noise in the FIR maps (see discussion in Sect. [^. The 
stacking maps are the results of median stacking the Herschel 
maps for all the five sources with proper centring and signal- 
to-noise ratio at each wavelength. In case of SPIRE maps, we 
stacked the cleaned maps, resulting from the first iteration of 
source extraction described above. The stacking map at lOOyum 
is shown in Fig. [^ and did not provide a detection. Correspond¬ 
ing stacking maps at 160yum, 250yum, 350yum, and 500yum, re¬ 
spectively, were similar and did not provide a detection either. 
We performed photometry on the stacked maps in the same way 
as described in Sect. [4. 1 [ for the individual maps. The resulting 
flux density uncertainties are summarised in the last column in 
TableEl 


4.3. Stacking in the near- and mid-infrared regime 

Since the six observed IFRS were also undetected in the SWIRE 
maps at 24 yum, we performed a stacking analysis at these wave¬ 
lengths to obtain mean properties of our IFRS population in the 
following analysis. 

At 3.6yum, all six IFRS were covered by SWIRE, but only 
four by the deeper SERVS survey. We median stacked the 3.6yum 
maps at the positions of the six sources. We carried out aperture 
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Table 3. Resulting point source flux density uncertainties for the six IFRS at the flve observed wavelengths. The last column shows the respective 
uncertainties in the stacked maps, resulting from a median stack of flve individual maps at each wavelength. 
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Fig. 3. Resulting stacking map obtained from a median stack of the flve 
IFRS at 100 pm. IFRS S749 was excluded because of the higher noise as 
discussed in Sect. The individual maps were centred on the position 
of the IFRS which is marked by a cross in the centre of the map. No 
detection has been found in the stacking map. 


photometry using an annulus of 1.9"in radius and a background 
ring between 3.8"and 7.6"and applied aperture corrections. No 
source was found in the stacked map. The uncertainty was ob¬ 
tained based on 225 randomly placed apertures and background 
rings of same size and iteratively removing those apertures re¬ 
sulting in flux densities exceeding 3cr. We found an uncertainty 
of (^ 3 . 6 fim = 0.397yuJy in the stacked map. 

At 24 pm, all six IFRS were covered by SWIRE. We median 
stacked the six individual maps and carried out aperture photom¬ 
etry using an annulus of 3.5"in radius and a background ring be¬ 
tween 6.0"and 8.0"and applying aperture corrections. No source 
was found in the stacked map either. We obtained the uncertainty 
in the same was as described for the stacked map at 3.6 pm and 
found it to be cr 2 Afim = 31.3 pfy. 


5. Modelling and analysis 

The Herschel observations between 100/im and 500/im of the 
six IFRS presented in this work did not provide a detection. A 
stacking analysis did not show an FIR counterpart below the de¬ 
tection limit either. However, the flux density upper limits mea¬ 
sured in Sect.j^put constraints on the SED of IFRS. 

The six IFRS observed with Herschel were selected to show 
high radio-to-IR flux density ratios (> 2000; see Sect.j^. Five of 
these IFRS are undetected in the near-IR regime while SERVS 
provides a 3.2(T couterpart at 3.6 pm for one IFRS as reported by 
|Maini et al.| ( |2013| ). This counterpart was unknown at the time 
when the observations with Herschel were carried out. Except 
for this single near-IR detection, these six IFRS are solely de¬ 
tected in the radio regime. Thus, no redshift is known for these 
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sources. By these characteristics, the six IFRS observed with 
Herschel are among the most extreme objects in the class of 


IFRS. 


All IFRS with known spectroscopic redshifts, presented by 
Collier et ar]p014| ) and |Herzog et al.| ( [T014| ), are faintly detected 


in the near-IR regime. Their redshifts, 1.8 < z < 3, are in agree¬ 
ment with the suggestion by |Norris et al.| ( |2QlT] ) that IFRS follow 
the correlation between 3.6yum flux density and redshift found 
for HzRGs by [Seymour et al.| ( |2007| ). Assuming that this cor¬ 
relation holds at lower IR flux densities, as it does for HzRGs, 
the IFRS observed with Herschel would be placed at redshifts 
of z > 4 or higher because of their 3.6 pm faintness. 


We emphasise that no redshift above z > 3 has been mea¬ 
sured for an IFRS, presumably because those are too faint for 
spectroscopic observations. Although we speculate that the IFRS 
discussed here are at z > 4, we will consider the broad redshift 
range 1 < z < 12. 


In the following, we tried to limit possible explanations for 
the phenomenon of IFRS, performing a broad-band SED mod¬ 
elling based on SED templates of known galaxies (Sect. |5.l| ). We 
started with the simplest approach and shifted these templates 
in a broad redshift range, testing them against the photometric 
data of IFRS (Sect. [5.1 H] ). Subsequently, we broadened the pa¬ 
rameter space by modifying our broad-band SED templates. We 
scale d the te mplates in luminosity (Sect. [5.1.^ , added extinction 
(Sect. [5.1 3] ), and finally modified them simultaneously in lumi¬ 
nosity and extinction (Sect. |5.l!4| ). In each approach, we tested 
their compatibility with the photometric constraints of IFRS . Fi¬ 
nally, we constrained the IR SED of IFRS based on decompos¬ 
ing their IR SED into an AGN and a starburst (SB) component 
(Sect. |5.2| ), both putatively contributing to the total emission of 
this peculiar class of objects. 


5.1. Broad-band SED modelling 

In order to model the broad-band SED of IFRS, we used pho¬ 
tometric data of our objects in all available wavelength regimes. 
The IFRS observed with Herschel provided 1.4 GHz fiux den¬ 
sities between 7 mJy and 26 mJy. In the modelling, we used a 
median 1.4 GHz fiux density of 15 mJy, but also discuss the out¬ 
come for higher and lower radio fiux densities. Furthermore, we 
used available constraints at 3.6pm and 2Apm. At 3.6pm, we 
used a median 3cr fiux density upper limit of \.\9p3y as dis¬ 
cussed in Sect. |4.3| In the mid-IR at 24 pm, we used a 3cr fiux 
densi ty up per limit of 94 p]y based on the stacking described in 
Sect. |4.3[ Furthermore, we made use of the far-IR fiux density 
upper limits at five wavelengths between lOOyum and 500 yum 
measured in this work. More precisely, we used 3cr fiux density 
upper limits based on the uncertainties in the median stacked 
maps presented in Sect. |4.2| and summarised in the last column 
in Table O 
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Based on photometric data and redshifts from the 
NASA/IPAC Extragalactic Database (NED), we built SED 
templates for different objects by connecting the data points 
and smoothing the template where appropriate. We used tem¬ 
plates of the spiderweb galaxy (an HzRG, also known as 
MRC 1138-262), the local radio galaxy Cygnus A (also known 
as 3C 405), the CSS source 3C 48, the RL quasar 3C 273, the lo¬ 
cal ULIRG Arp 220, the Seyfert galaxy Mrk 231, the local SB 
galaxy M82, an RL hyper-luminous infrared galaxy (HyLIRG; 
IRAS E15307-f3252), IRAS E00183-7111 (referred to as 
E00183; this object is a ULIRG, showing contribution of an 
RL AGN combined with significant SB activity), and the qui¬ 
escent elliptical brightest cluster galaxy NGC 1316. We added 
photometric IR data between 100yum and 870 yum from |Seymour| 
|et al.| ( |2012| ) to the template of the spiderweb galaxy since the 
IR coverage of this template is poor in NED and this wavelength 
regime is crucial for our analysis. All photometric data points 
from [Seymour et al^ are at least 4cr detections. 


5.1.1. Shifting broad-band SED templates to various 
redshifts 


As a simplest approach, we tested whether any SED template of 
a known galaxy is in agreement with the available photometric 
data of lERS when these templates were shifted to different red¬ 
shifts, keeping constant rest-frame luminosity and scaling flux 
densities using a flat ACDM cosmology (Qa = 0.7, Qm = 0.3, 
Hq = 70kms“^ Mpc“^). The resulting plots for redshifts of 1, 4, 
8 , and 10, respectively, are shown in Eig.[^ 

In the redshift regime 1 < z < 5, no template fulfils all avail¬ 
able photometric constraints, i.e. none of these objects would 
produce the observational characteristics of lERS if placed at 
redshifts z < 5. However, for redshifts in the range 5 < z < 8.5, 
Cygnus A provided an appropriate template to fulfil the pho¬ 
tometric constraints of lERS. The spiderweb galaxy SED was 
in agreement with all available photometric data of lERS when 
shifted to z > 10.5. At lower redshifts, the radio and the near-, 
mid-, and far-IR flux densities of the spiderweb galaxy exceeded 
the measured fluxes of lERS. We also tried templates of other 
HzRGs (e.g. 3C470, 4C 23.56) and found the same qualitative 
result. 

The CSS source 3C48 and the RL quasar 3C273 matched 
the measured radio flux density of lERS if shifted to the redshift 
regimes 6 < z < 9 and 8 < z < 12, respectively. However, at 
these redshifts, their 3.6yum flux densities exceeded the related 
measured flux density upper limits of lERS up to a factor of ten, 
ruling out these templates. All other templates were found to be 
in strong disagreement with the lERS data at all redshifts be¬ 
tween 1 and 12. Either the radio flux densities of the templates 
exceeded the measured flux density of lERS or the templates dis¬ 
agreed with the near-, mid-, and far-IR flux density upper limits 
of lERS by several orders of magnitude. 

We also estimated a score for each SED template at each red¬ 
shift. Eor this purpose, we derived the score 5' based on a modi¬ 
fied Li norm (e.g. Horn & Johnson|19'85] ) given by 


10 


- 2 logio - 2 H[mj-3crj] logjo 


( 1 ) 


The first sum in the exponent runs over the wavelengths i at 
which our sample of lERS is detected, where y denotes the 
flux density of lERS and m the flux density of the respective 
SED model. The second sum accounts for deviations of the SED 
model from the flux density upper limits of lERS and sums over 
the wavelengths j at which the lERS are undetected with an 


rms (T, using a Heaviside step function H. In order to avoid an 
overweighting of the EIR regime when deriving a score due to 
the five independent flux density measurements in this regime, 
we only used one out of the five Herschel wavelengths. Thus, 
the wavelength resulting in the highest addend \ logiQ(3crj/mj)\ 
was used for this purpose. 

This implies that a score of 1 is obtained by an SED which is 
in agreement with all eight photometric constraints of the lERS 
considered in this work. We emphasise that a score of 1 does not 
imply that the SED of lERS necessarily follows the respective 
template. However, the score provides useful insights how well 
a template agrees with the data of lERS at a respective redshift. 
The upper left subplot in Eig.|^ shows the score as a function of 
redshift for different templates. 

We found only the spiderweb galaxy and Cygnus A are con¬ 
sistent with the photometric data of lERS when shifted to appro¬ 
priate redshifts. If lERS have the SED of a known galaxy, lERS 
have to be at redshifts z > 5. This is consistent with the sug- 
gestion by|Norris et al.| ( [201l] ), |Collier et al.| ( |2Q14| ), and IHerzo'g 
|et al.| ( |2014| ) that the IFRS observed with Herschel are at red¬ 
shifts z > 4 because of their faintness at 3.6yum as discussed 
above. 


5.1.2. Modifying broad-band SED templates in redshift and 
luminosity 


In the broad-band SED modelling presented in Sect. |5.1.1| red¬ 
shift was the only free parameter. We found that only HzRGs 
(e.g. the spiderweb galaxy) and Cygnus A provide appropriate 
templates to explain the emission characteristics of IFRS, and 
only if shifted to redshifts z > 5. In the following, however, we 
aimed at testing whether the SED characteristics of IFRS can be 
reproduced by objects at lower redshifts. For this purpose, we 
broadened the parameter space and used templates modified by 
varying their luminosity. 

In this approach, we included a wavelength-independent lu¬ 
minosity scaling factor as an additional parameter. This allowed 
the templates to be scaled to match the mean radio flux density 
of the IFRS observed with Herschel (5 'i.4ghz ~ 15mJy), since 
these sources are solely detected in the radio regime apart from 
the weak 3.6yum counterpart of S427. 

The upper right subplot in Fig. shows the score versus red¬ 
shift for the different templates when they were scaled in lu¬ 
minosity to match the median radio flux density of the IFRS 
discussed in this work. In this approach, we found the spider¬ 
web galaxy is in agreement with the available photometric con¬ 
straints at redshifts in the range 4 < z < 6.5 when scaled in 
luminosity by factors between 0.08 and 0.4. At higher redshifts, 
the scaled template SED exceeded the 24 yum flux density upper 
limits and therefore resulted in a score below I. The SED tem¬ 
plate of Cygnus A was also found to be in agreement with the 
data at z ~ 1 and in the redshift range 4 < z < 10 when scaled in 
luminosity by factors between 0.007 and 2. In the redshift range 
1.5 < z < 4, the scaled template disagreed with the 3.6yum flux 
density upper limit. All other templates were in disagreement 
with at least one upper limit at any redshift. 

When scaling the templates to match the faint end of the 
observed radio flux density range of IFRS of 7 mJy, the spider¬ 
web galaxy template was in agreement with all constraints when 
shifted to the redshift range z > 3.5, and Cygnus A if shifted to 
the range z < 2 or 3 < z < 10.5. None of the other templates 
agreed with all photometric constraints. In contrast, if the tem¬ 
plates were scaled to the radio flux density of the radio-brightest 
IFRS in our sample (26 mJy), only the template of Cygnus A was 
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Fig. 4. Broad-band SED modelling for IFRS, using the far-IR flux density upper limits measured in this work (red arrows), SERVS/SWIRE 
near- and mid-IR flux density upper limits (black arrows), and the detections at 1.4 GHz (black bar). The template SEDs were shifted to redshifts 
between z = 1 and z = 12. Shown are the results at z = 1 (upper left subplot), z = 4 (upper right), z = 8 (lower left), and z = 10 (lower right), 
respectively. The legend shown in the upper left subplot is valid for all four subplots. At redshifts z < 5, no template is in agreement with the 
photometric constraints of IFRS. However, if shifted to the redshift regime 5 < z < 8.5, Cygnus A (orange solid line) is in agreement with all 
constraints. The spiderweb galaxy (black dashed line) provides an appropriate template if shifted to very high redshifts z > 10.5. The CSS source 
3C 48 (blue dotted line) and the RL quasar 3C 273 (magenta dashed-dotted line) fulfil the constraint in the radio regime at redshifts 6 < z < 9 and 
8 < z < 12, respectively, but disagree with the 3.6 fim flux density upper limits at these redshifts. 


in agreement with the IR flux density upper limits (4.5 < z < 10). 
However, at z ~ 4, the template of the spiderweb galaxy only 
slighly violated the 24 fim flux density upper limit. 


5.1.3. Modifying broad-band SED templates in redshift and 
extinction 


We found in Sect. |5.1.1| that the redshifted SED templates of 
3C 48, and 3C 273 were in agreement with the radio flux density 
range and the FIR flux density upper limits of IFRS, however ex¬ 
ceeding the near- or mid-IR flux density upper limits. Therefore, 
here, we tested whether some of our templates could be squeezed 
below these IR flux density upper limits when additional extinc¬ 
tion was added to the templates. 

For this purpose, we added extinction in the rest-frame opti¬ 
cal and near-IR regime of the template SEDs if required. We used 
the |Calzetti et al.|p000| ) reddening law which reduces emission 
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at rest-frame wavelengths between 0.12 jjm and 2.2 jjm. We lim¬ 
ited the additional extinction to an arbitrary number of 500 mag. 
In our modelling, we applied energy conservation by calculating 
the power hidden by the additional amount of dust in the rest- 
frame optical and near-IR regime and re-radiating this power at 
mid- and far-IR wavelengths. We implemented this power con¬ 
servation in our modelling by adding a dust emission component 
of the same power as hidden by the additional extinction in the 
optical and near-IR regime. 

The emission of dust is usually described by a modified 
black-body spectrum, i.e. Sy oc By(T)6y, where By(T) is the 
black-body emission given by the temperature T and 6y is the 
emissivity function at a frequency v. We followed simple ap¬ 
proaches, using only one dust temperature T = 70K and as¬ 
suming e oc yP with p = 1.5. The emissivity spectral index p 
depends on dust grain properties as size and composition and 
usually ranges between 1 and 2. We note that the choice of this 
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Fig. 5. Score for different SED templates as a function of redshift, measured in our broad-band SED modelling. A score of 1 implies that the SED 
is in agreement with all photometric constraints of the lERS discussed in this work. Shown is the score resulting from modelling without scaling 
and without additional extinction (upper left), with scaling and without additional extinction (upper right), without scaling and with additional 
extinction (lower left), and with scaling and with additional extinction (lower right subplot), respectively. In both plots on the right hand side, the 
templates were scaled to a 1.4 GHz flux density of 15 mJy. The legend shown in the upper left subplot is valid for all four subplots. 


spectral index in the given range did not qualitatively change our 
results; neither did changing the dust peak temperature between 


30 K and 100 K. 


The plot of the scores resulting from this modelling 
approach—adding extinction to the redshifted SED templates— 
is shown in the lower left subplot in Fig. It should be noted 
that no scaling in luminosity as presented in Sect. |5.1.2] had been 
applied in this approach. We found the spiderweb galaxy tem¬ 
plate in agreement with the data at redshifts z>9.5. Atz~ 10, 
5.6 mag of additional extinction brought this SED template in 
agreement with the data. At higher redshifts, no additional ex¬ 
tinction was required. Cygnus A was in agreement with all flux 
density upper limits in the redshift range 5 < z < 9 without 
adding extinction as shown in Sect. |5.1.T] 3C273 matched the 
data at redshifts z > 7.5 if between 8.6 mag and 1.4 mag of ex¬ 
tinction were added. 3C48 was in agreement with all IR flux 
density upper limits at z > 7 if between 1.4 mag and 0.2 mag of 
extinction were added. However, the radio flux density of 3C 48 
deceeded the minimum flux density of 7 mJy at high redshifts 


z > 9.5. All other SED templates were found to be in disagree¬ 
ment with the data at any redshift. 

Even enormous and unphysically high amounts of extinc¬ 
tion of 500 mag were insufficient to squeeze the near-IR emis¬ 
sion of some templates below the flux density upper limit at the 
observed wavelength of 3.6)um. In that case, we did not add ex¬ 
tinction to the templates, resulting in the appropriately low score 
at the respective redshift. This happened for all radio-quiet tem¬ 
plates like Arp220, Mrk231, or M82 at redshifts z < 5. How¬ 
ever, at higher redshifts, when additional extinction was able to 
squeeze the templates below the IR flux density upper limits, the 
templates were by far too faint in the radio regime to explain the 
emission features of IFRS. 


5.1.4. Modifying broad-band SED templates in redshift, 
luminosity, and extinction 

Finally, in the last step of our broad-band SED modelling, we 
combined the approaches presented in Sections [5.1.2| and [5. 1.3 
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Here, we modified the SED templates by shifting them in the 
redshift range 1 < z < 12, scaling them in luminosity to 
match the observed radio fiux density, and adding extinction— 
if required—to squeeze the SED templates below the near- and 
mid-IR fiux density upper limits. IR re-radiation was imple¬ 
mented as described in Sect. 15.1.31 

The resulting plot of the score as a function of redshift is 
shown in the lower right subplot in Eig. We note that each 
template at each redshift had been modified by an individual 
scaling factor and—if necessary—an individual amount of ex¬ 
tinction had been added. 

Similar to the approaches presented before, we only found 
the spiderweb galaxy, Cygnus A, the CSS source 3C48, and 
the RE quasar 3C273 to provide suitable templates to explain 
the observational data of lERS. The spiderweb galaxy was in 
agreement with the data at all redshifts except for the range 
6.5 < z < 9.5. In this redshift range, the template exceeded the 
24 pm fiux density upper limit despite adding additional extinc¬ 
tion. The template of Cygnus A matched all constraints at z < 10 
and 3C 273 at redshifts z > 2.5. 3C 48 was in agreement with the 
data at any redshift except for very low redshifts z < 1.5 and 
except for the range 2.5 < z < 7. 

In the respective redshift ranges, an appropriate scaling fac¬ 
tor and an appropriate amount of extinction was found for the 
spiderweb galaxy, Cygnus A, 3C48, and 3C273 to match all 
photometric constraints known for these lERS. Eor example at 
z = 4, the spiderweb galaxy and Cygnus A needed to be scaled 
down in luminosity by a factor of 10 and 3, respectively, to match 
the radio fiux density of 15mJy. No extinction was needed in 
these cases. In contrast, 3C273 required 3.2mag of extinction 
when scaled down by a factor of 5 to match the radio data point, 
while 3C 48 could not be modified at this redshift to match all 
constraints. 

Matching the radio fiux densities of the radio-brightest lERS 
in our sample (26mJy at 1.4 GHz), again in case of a shift to 
z = 4, an additional extinction of 1.3 mag was required in case 
of Cygnus A, whereas no reasonable amount of extinction could 
be added to the templates of the spiderweb galaxy, 3C48, and 
3C 273, to match the near- and mid-IR fiux density upper limits. 
Considering the radio-brightest lERS in our sample, at z < 4.5, 
only the Cygnus A SED template could be brought in agreement 
with the constraints of these most extreme—i.e. with the highest 
radio-to-IR fiux density ratios—lERS, if up to 4.1 mag of extinc¬ 
tion were added. At higher redshifts z > 10, 3C48 and the spi¬ 
derweb galaxy provided appropriate templates, too. 3C273 was 
in agreement with the data at 4.5 < z < 5.5 and 7.5 < z < 12. In 
contrast, if the templates were scaled to a 1.4 GHz fiux density 
of 7mJy, the spiderweb galaxy, Cygnus A, 3C273, and 3C48 
could be modified to be in agreement with all photometric con¬ 
straints at any redshift, except for Cygnus A at very high red¬ 
shifts z > 10. 

All other templates used in this broad-band SED modelling 
disagreed with the observational data when shifted to the red¬ 
shift range 1 < z < 12, scaled in luminosity, and extinction 
was added. These templates were clearly ruled out by their dis¬ 
agreement with the fiux density upper limits of lERS in the EIR 
regime measured in this work. Particularly, the templates of an 
RE HyLIRG or E00183, representing common classes of lumi¬ 
nous objects at z > 2, were found to be inappropriate to repro¬ 
duce the SED constraints of lERS at any redshift. 
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5.1.5. Summary of broad-band SED modelling 


We summarise our findings from the broad-band SED mod¬ 
elling as follows. When using templates of existing galaxies 
(Sect. |5.1.1| ), the SED characteristics of lERS can only be ex¬ 
plained by the SEDs of Cygnus A and HzRGs, however only at 
high redshifts 5 < z < 8.5 and z > 10.5, respectively. If lERS are 
at lower redshifts (z < 5), their SED constraints can be fulfilled 
by fai nter versions of an HzRG or Cygnus A (Sects. [5X^ and 
|5.1.4| ), by a dust-obscured CSS source, or by a dust-obscured 
RL quasar (Sects. [5X^ and [5X4| ). 

If lERS are linked to HzRGs as suggested by [Middelber^ 
et al.| (|2Q11|), INorris et al.| ( |2Q11| ), |Collier et al.| ( |2014| ), and 


Herzog et al.| (2014), there are two different suggested op¬ 


tions: (a) lERS are very similar to HzRGs—i.e. with similar 
luminosities—however at higher redshifts (z > 5), or (b) lERS 
are fainter siblings of HzRGs at similar redshifts (1 < z < 5). 
Based on the template of the spiderweb galaxy, we found that 
both options are consistent with the EIR fiux density upper lim¬ 
its of lERS measured in this work and the related broad-band 
SED modelling. If lERS have the same intrinsic properties as 
HzRGs, represented by no scaling and no adding of extinction 
in our modelling, they have to be located at very high redshifts 
z > 10.5 (Eig. upper left subplot). Alternatively, if lERS are 
fainter siblings of HzRGs, i.e. they are scaled down in luminos¬ 
ity and potentially more dusty, they could be located at lower 
redshifts z > 1 (Eig. lower right subplot). It should be noted 
that extinctions between 14 mag and 5 mag had to be added to 
the scaled spiderweb galaxy SED template at redshifts z < 2 to 
match the IR fiux density upper limits. 

Two of the templates, 3C 48 and 3C 273, were found to be in 
agreement with the data of lERS only if additional extinction in 
the order of up to several magnitudes was added to the templates 
as discussed in Sect. |5.1.4| Particularly, these modified templates 
did not exceed the EIR fiux density upper limits of lERS mea¬ 
sured in this work. This implies that a significant dust obscu¬ 
ration in lERS, explaining the optical and near-IR faintness of 
these objects, cannot be ruled out. In fact, if the lERS observed 
here are similar to 3C 48 or 3C 273, these sources have to be dust 
obscured to produce the photometric constraints of lERS. 

The results of this broad-band SED modelling are generally 
in agreement with the findings by [Herzog et al.| ( 20 14| ). [Herzog | 
[et al.[ measured spectroscopic redshifts of three lERS detected in 
the optical and near- and mid-IR regime and did the first redshift- 
based SED modelling for the class of lERS. They find the mod¬ 
ified templates of 3C 48 and 3C 273 in agreement with all avail¬ 
able photometric data of these three lERS, whereas the templates 
of star forming galaxies and Seyfert galaxies were clearly ruled 
out. We also agree with the finding that only RL obj ects can 
explain the photometric constraints of lERS. However, [Herzog 


[et al.[ find that no additional extinction is required for 3C 48 and 


3C 273 in their SED modelling. Here, in contrast, we found that 
lERS might be dust obscured RL quasars with spectra similar to 
3C 48 or 3C 273. However, [Herzog et~ar] studied the IR-brightest 
lERS, whereas the IR-faintest lERS are studied i n the current 
work, providing a valid reason for the discrepancy. [Huynh et~ar 


( [2010| ) modelled the SED of individual lERS and find that the 
scaled SED template of 3C273 is in agreement with their data 
if around one magnitude of extinction was added. This in agree¬ 
ment with our modelling, considering that the lERS from [Huynh[ 
et al.[ are more than ten times fainter in the radio compared to 


the lERS discussed in the present work. We found that lERS 
have to be dust obscured, too, if lERS are similar to HzRGs at 
z < 3.5. However, if the spectral shape of lERS is related to that 
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of HzRGs, no extinction is required to explain the SED charac¬ 
teristics of IFRS if these templates are shifted to higher redshifts. 

Our broad-band SED modelling showed that the SED tem¬ 
plates of the spiderweb galaxy and of Cygnus A can be shifted 
to I < z < 2 and modified to agree with the photometric ob¬ 
servations of IFRS. However, in contrast, all 21 spectroscopic 
redshifts of IFRS are in the range 1.8 < z < 3, i.e. no IFRS 
has been found at lower redshifts z ~ 1. There are two poten¬ 
tial reasons for this discrepancy, (a) Our modelling—by shifting 
the templates in redshift, scaling them in luminosity, and adding 
extinction—is not physically realistic for galaxies at z ~ 1. At 
these low redshifts, all appropriate templates needed to be scaled 
down in luminosity up to factors of 500. Presumably, the black 
hole mass needs to be scaled down by a similar factor. How¬ 
ever, RE AGN with lower-mass black holes are very rare (Rees 
et al., submitted). If the galaxy loses its radio excess by scaling 
down the black hole mass, the characteristically high radio-to-IR 
flux density ratio of IFRS would not be reached and the galaxy 
would not be considered as an IFRS. Furthermore, in case of the 
spiderweb galaxy template, the required high amounts of addi¬ 
tional extinction of around 10 mag at low redshifts are very rare 
and decrease the possibility to observe such extreme objects, (b) 
The IFRS discussed in this work are intrinsically different to the 
IFRS with known spectroscopic redshifts. Cygnus A provided 
the only template which can reproduce the characteristics of 
IFRS at z ~ 1 without additional dust obscuration. However, it is 
known that Cygnus A contains a hidden quasar ( [Antonucci et al.| 
I994| ), resulting in narrow emission lines, whereas IFRS were 


found to show broad emission lines in their optical spectra ( |Col-| 
lier et al.|20I4t|Herzog et al.|20I4| ). If the IFRS discussed in this 
work are indeed fainter versions of Cygnus A at redshifts z ~ 1, 
these IFRS would have different properties than the IR-detected 
IFRS presented by |CoIlier et al.| and [Herzog et al.| They would 
form a separate subclass with different redshift and emission line 
properties. However, so far, no evidence has been found that the 
population of IFRS might be devided into two sub-classes. 

We note that the modified templates of Cygnus A, 3C 48, and 
3C273—found to match all photometric constraints of IFRS— 
can be considered as HzRGs in our modelling at z > 1. At these 
redshifts, they fulfil both selection criteria of HzRGs: z > 1 
and L 3 GHZ > 10^^ WHz“^ Although we scaled down the tem¬ 
plates in luminosity at lower redshifts, they were radio-luminous 
enough to be considered as HzRGs. We can therefore conclude 
that only HzRGs can explain the photometric characteristics of 
IFRS. 

Our broad-band SED modelling showed several templates 
that contain a significant SB contribution in disagreement with 
the available photometric data (Arp220, M82, RF HyFIRG, 
F00I83). In contrast, the highly star forming spiderweb galaxy, 
fulfilled the photometric data constraints of IFRS. Therefore, in 
the following, we analysed a potential star forming contribution 
to the SED of IFRS, instead or in addition to an RF AGN. 


5.2. IR SED modelling 

One of the crucial questions in revealing the nature of IFRS con¬ 
sists in unmasking the star forming activity and its contribution 
to the energy output in contrast to the activity of the galactic nu¬ 
cleus. In the IR regime, the emission of galaxies is mainly given 
by these two components, i.e. AGN and SB. The AGN emission 
comes from dust, reprocessing far-UV through optical light, and 
peaks at around lOyum, corresponding to a dust temperature of 
300 K. In contrast, the stellar component is strongest in the FIR 
and arises from stellar emission reprocessed to the FIR regime 


by dust grains with a maximum at a wavelength around lOOyum, 
corresponding to 30 K. 


5.2.1. IR SED fitting based on IR flux density upper limits 

In order to decompose the IR emission of HzRGs, [Prouart et al. 
( |20I4| ) fit the IR SED of HzRGs constrained by Herschel and 
ancillary IR data. They used an IR SED model based on the as¬ 
sumption that the IR emission comes from dust heated by star 
forming activity and an AGN. They built template SEDs, adding 
an empirical AGN template and one SB template out of a set of 
empirical SB templates. 

Here, we followed this approach, aiming at setting upper lim¬ 
its on the IR SED of IFRS based on the available FIR flux density 
upper limits measured in Sect. We used templates derived by 
|Mullaney et al.| ( [MTT] ), cove ring the rest-fram e wavelength range 
between 6 yum and 1090 yum. |Mullaney et al.l empirically built five 
different SB templates, differing in their peaking temperature as 
well as in the strength of the emission of the polycyclic aromatic 
hydrocarbon (PAH) molecules. These templates cover the entire 
range of host galaxies (see [Mullaney et al.||2QTT] for a detailed 
description of the SEDs). Based on these SB templates, they de¬ 
rived one AGN template as the residual SED after removing the 
SB contribution. 

We used these templates and built the total IR SED by mul¬ 
tiplying the AGN template by a wavelength-independent lumi¬ 
nosity scaling factor and adding one SB template which was also 
multiplied by a wavelength-independent luminosity scaling fac¬ 
tor. We set an upper limit on the IR emission by modelling the 
SED in order to minimise the deviation between the total SED, 
i.e. the sum of SB and AGN template, and the available observed 
3(T fir flux density upper limits from the stacked maps (see last 
column in Table [^. It was required that the available flux den¬ 
sity upper limits were fulfilled in the modelling. In this proce¬ 
dure, the most appropriate SB template out of the five available 
templates was determined based on the lowest deviation. At red¬ 
shifts z < 3, the observed 24/im emission fell into the templates’ 
wavelength coverage. At these redshifts, we use d the 3(T flux 
density upper limit at 24yum—measured in Sect. |4.3| from the 
stacked map—in addition to the FIR flux density upper limits 
measured in the present work. 

By this approach, we determined exactly that total IR SED— 
consisting of the AGN template and one SB template—which 
has maximum IR flux but is still in agreement with all avail¬ 
able flux density upper limits. Figure shows the resulting 
IR SED modelling for z = 4. We emphasise that our modelling 
represents an upper limit on the total IR emission of IFRS. 

This approach implies that the true AGN contribution to the 
emission of the IFRS could be higher than the contribution cal¬ 
culated from the maximum total model described above, though 
accompanied by a lower SB contribution, and the other way 
around. However, the difference to the highest possible contri¬ 
bution of each component is generally rather low. 

By this modelling, aiming at setting upper limits on the 
IR emission and referred to as “maximum model”, we decom¬ 
posed the maximum IR SED into an AGN component and an 
SB component. Based on this decomposition, we were able to 
derive IR luminosities of both components by integrating the 
flux density in the rest-frame wavelength range between 8 jum 
and lOOOyum. Figureshows the IR luminosities versus redshift 
for both components and the total IR luminosity of the maximum 
model. Furthermore, the figure shows the absolute upper limit 
for both the SB and the AGN IR luminosity, assuming that the 
other IR emitting process does not contribute at each time. For 
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Fig. 6. IR SED modelling for IFRS at redshift z = 4, using the far- 
IR flux density upper limits (red arrows) measured in this work. The 
total (black solid line) IR emission is composed of an AGN compo¬ 
nent (red dotted line) and a SB component (blue dashed line). The max¬ 
imum model represents an upper limit on the IR emission of IFRS and 
is deflned by the highest IR flux which is in agreement with all available 
flux density upper limits. 



Fig. 7. Infrared luminosity versus redshift for the SB (blue) and 
AGN (red) component as well as for the total IR luminosity (black). The 
dashed lines show the data of the maximum IR SED model, combining 
the AGN and the most appropriate SB template. The solid lines show 
the IR luminosity upper limits for both components, assuming that the 
other component does not contribute to the IR emission at each time. 


both components, the luminosities of the maximum model are 
generally only slightly lower than the absolute upper limit lumi¬ 
nosities, except for the SB component at very high redshifts. We 
note that the jump in the maximum IR luminosity of the AGN 
and of the SB component at z = 3 is related to the flux density 
upper limit at 24pm which only coincides with the template’s 
wavelength range at redshifts z < 3. 

We emphasise that almost any decomposition of the IR SED 
into SB and AGN component is possible for IFRS since this 
wavelength regime is only constrained by upper limits. There- 
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fore, we are only able to set upper limits on the luminosity. Par¬ 
ticularly, we note that the maximum model shown in Figures 
and is not more likely than any other combination of SB and 
AGN components which is in agreement with the flux density 
upper limits. However, the maximum model represents exactly 
that combination of SB and AGN components which results in 
the highest total flux but is still in agreement with all flux den¬ 
sity upper limits. Therefore, this maximum model sets the upper 
limit on the total IR luminosity. 

We found that the FIR flux density upper limits measured 
in this work constrain IFRS to have total IR luminosities of 
< 10^^ Lq at redshifts z < 1.5. At z < 2.5, IFRS can have a max¬ 
imum IR luminosity between 10^^ Lq and 10^^ Lq. Such sources 
are labelled as UFIRGs. At redshifts z < 6, the IR luminosity of 
IFRS can be at most 10^^ L© which allows IFRS to be HyFIRGs. 
Drouart et al. | ( |20 1 4\ And absolute numbers and upper limits 


for the total IR luminosity of HzRGs in the range of a few to a 
few tens of 10^^ L©. Our redshift-dependent IR luminosity upper 
limits of IFRS are in agreement with these numbers of HzRGs. 

We note that the maximum model shown in Fig. |7] (dashed 
lines) is unphysical for redshifts z < 2.5. For these redshifts, 
the AGN contribution to the maximum model is very low, cor¬ 
responding to a very low AGN activity or even an absent AGN. 
Instead, according to that model, the IR emission is mainly or 
even completely produced by SB activity. In case of a purely 
star forming galaxy, a direct connection between the radio and 
the IR emission was found (i.e. |Yun et al.pOOT] ), known as radio- 
IR correlation. Using the outcome of the maximum model that 
the emission of IFRS is dominated by the SB component for 
redshifts z < 2.5, we could estimate the radioluminosity from 
the IR luminosity. Using our FIR luminosity upper limits, we 
found that this radio luminosity is at least two orders of magni¬ 
tude lower than the radio luminosity calculated from the mea¬ 
sured 1.4 GHz flux density of IFRS. This discrepancy clearly 
showed that the existence of an AGN in the IFRS is essential 
to explain their radio flux densities. Modelling their IR SED by 
only an SB template is inappropriate. The existence of an AGN 
in IFRS is also in agreement with VFBI obse rvations by [Norris 


|et al.| (|2007|), |Middelberg et al.| (|2008b|), and |Herzog, A. et al. 


5.2.2. Including the radio data point into the IR SED 
modelling 

The IR SED modelling presented above, aiming at setting upper 
limits on the total IR flux of IFRS, led to unphysical results be¬ 
cause it neglected the measured radio flux density of IFRS. In 
this section, we attempted to address this issue by expanding our 
modelling to the radio regime based on an RF AGN template. 

In Sect. im we showed that all SEDs in agreement with the 
photometric characteristics of IFRS were HzRGs and that the 
spiderweb galaxy provided an appropriate SED template in the 
broad-band SED modelling. It is known that both star forming 
activity and an AGN contribute to the emission of this power¬ 
ful galaxy. Therefore, in the following, we used the broad-band 
SED of the spiderweb galaxy, built from NED and the additional 
IR data from [Seymour et aL]p012| ) as described in Sect. [5.1 1 as 
a basic template, modifled by a wavelength-independent lumi¬ 
nosity scaling factor. Beyond that, to account for an SB com¬ 
ponent independent of the component already present in the 
spiderweb galaxy template, we added the most app ropriate SB 
template from the sample by jMullaney et al.jpOllj ). Similar to 
our approach in Sect. [5.2. 1| we modified this SB template by 
a wavelength-independent luminosity scaling factor. However, 
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Fig. 8. IR SED modelling for IFRS at redshift z = 4, using the far- 
IR flux density upper limits (red arrows) measured in this work, the 
SERVS/SWIRE flux density upper limits (black arrows), and the mea¬ 
sured radio flux density of 15mJy at 1.4 GHz. The total (solid line) 
emission is composed of the scaled spiderweb galaxy template (dotted 
line) and one scaled SB component (dashed line). Note that the limited 
wavelength coverage of the SB template caused the leaps at SOjum and 
5450/rm. The fit represents an upper limit on the IR emission of IFRS 
and is deflned by matching the radio flux density of 15 mJy at 1.4 GHz 
and providing the lowest deviation between all available flux density 
upper limits and the total SED. Note that the dash at 21 cm (flux density 
of ^ 20 jLiJy) is the contribution of the additional SB component to the 
1.4 GHz flux density. 


here, we also used the radio-IR correlation ( |Yun et al.|2001| ) to 
estimate the radio luminosity of the star formation component, 
in addition to that from the RL AGN. 

R ecapit ulating, similar to the approach presented in 
Sect. |5.2.1| we composed our total SED of two different com¬ 
ponents. The spiderweb galaxy contributed SB and AGN emis¬ 
sion, while the additional SB component contributed additional 
SB emission. Based on the two independent scaling factors, both 
the AGN contribution and the SB contribution could be varied 
independently. We note that the scaling factor of the additional 
SB component could also be negative, i.e. representing a lower 
star forming activity than that in the spiderweb galaxy. 

Similar to our approach in Sect. |5.2.f| we modelled the max¬ 
imum IR SED by maximising the flux of the total template—i.e. 
composite from the spiderweb galaxy template and one addi¬ 
tional SB template—however, requesting the template to be in 
agreement with all available flux density upper limits. Eurther- 
more, here, we required the total template to match the median 
measured 1.4 GHz flux density of the observed lERS of 15 mJy. 
Eigurej^ shows the resulting maximum model for z = 4. 

We found that the contribution of the additional SB com¬ 
ponent to the radio emission, calculated based on the radio-IR 
correlation, is negligible compared to the radio emission of the 
spiderweb galaxy template. In case of z = 4, this SB contribution 
is around 20 juJy at 1.4 GHz which is more than two orders of 
magnitude lower than the contribution of the spiderweb galaxy 
template at this frequency. 

Based on the results of this modelling, we calculated the 
maximum total IR luminosity by integrating the flux density of 
the total SED template between 8 jum and lOOOyum in the rest 
frame. Eurthermore, we decomposed the maximum IR luminos- 



Fig. 9. Maximum infrared luminosity versus redshift for the AGN (red 
dotted line) component, the SB (blue dashed line) component, and for 
the total IR luminosity (black solid line). The underlying model (see 
Fig. for an example) matched the measured radio flux density of 
15mJy at 1.4 GHz. We note that the total IR luminosity and the SB 
IR luminosity are upper limits. 


ity in the SB contribution and the AGN contribution. |Seymour| 
et al. ( 2Q12| ) studied the IR properties of the spiderweb galaxy 


and And that 59% of its IR luminosity is contributed from the 
AGN and 41% from the star forming activity. We used these 
fractional numbers and decomposed the contribution of the spi¬ 
derweb galaxy template correspondingly. The maximum total IR 
luminosity and its decomposition in AGN and SB contribution 
as a function of redshift are shown in Fig. We found that, 
based on this modelling, IFRS are constrained to IR luminosi¬ 
ties below IO^^Lq at redshifts z < 2. At z < 5, IFRS can have a 
maximum luminosity of lO^^L©. At higher redshifts, we cannot 
exclude IFRS to be HyLIRGs which are deflned by IR luminosi¬ 
ties above IO^^Lq. 

We note that the luminosity upper limits obtained in this ap¬ 
proach are up to one order of magnitude lower than those in¬ 
ferred in Sect. |5.2.r] and presented in Fig. [7] This shows that it 
is crucial to include the measured radio flux density in the IR 
SED modelling. A modelling just based on the IR flux density 
upper limits is physically inappropriate and results in excessive 
IR luminosities. 


The numbers of the maximum total IR luminosity and the 
SB contribution are upper limits since they were obtained by a 
model aiming at the lowest deviation to the available IR flux 
density upper limits. In contrast, the IR luminosities of the AGN 
component can be considered to be rough estimates because of 
their direct link to the measured radio flux density. However, 
we point out that the IR emission of AGNs at high redshifts is 
known to spread by more than one order of magnitude ( |Drouart 
|et al.|20I^ . It should be emphasised that we assumed IFRS to be 
similar to HzRGs in this approach. Although several indications 
for this similarity have been found so far, it is not proven yet. 
While the luminosities shown in Fig. can be considered to be 
more realistic numbers, the IR luminosity upper limits presented 
in Fig.[7]are absolute upper limits which should not be exceeded, 
independent of the nature of IFRS. 

In our broad-band SED modelling described in Sect. mi we 
found that the local radio galaxy Cygnus A, the local CSS source 


Article number, page 13 of 


15 









































































ASlA proofs: manuscript no. document 


3C48 and the RL quasar 3C273 provide appropriate templates 
to explain the SED characteristics of IFRS, too. As mentioned 
above, these modified templates shifted to z > 1 fulfil the se¬ 
lection criteria of HzRGs. Therefore, we also used the SEDs of 
these galaxies as basic templates, i.e. containing an AGN and 
potentially a certain contribution from star forming activity, in¬ 
stead of the spiderweb galaxy. We performed the modelling in 
the same way as described above, i.e. using the scaled SED of 
Cygnus A, 3C 48 and 3C 273, respectively, as basic template and 
adding the most appropriate and scaled SB template, matching 
all available flux density upper limits and the measured 1.4 GHz 
flux density. However, in case of Cygnus A, 3C48 and 3C273, 
we did not find appropriate numbers in the literature to convert 
the IR emission of these objects into an AGN and an SB contri- 
butio n as we did for the sp iderweb galaxy based on the results 
from [Seymour et aL] ( |2Q12[ ). Therefore, we could only compare 
the total IR luminosity upper limits. We found that these num¬ 
bers derived based on Cygnus A, 3C 48 and 3C 273 as basic tem¬ 
plates are very similar to the numbers based on the spiderweb 
galaxy template shown in Fig. differing by not more than a 
factor of two. This finding provides evidence that the total IR 
luminosity upper limits presented in Fig. are indeed general 
upper limits, independent on the assumed nature of IFRS. 


5.2.3. Estimates of star formation and black hole accretion 
rate upper limits 

Based on the maximum IR luminosities measured in Sect. |5.2.2) 
we derived upper limit estimates of the star formation rate (SFR), 
causing the IR emission given by the SB component, and esti¬ 
mates of the black hole accretion rate, generating the AGN con¬ 
tribution to the IR flux. Since these numbers can only be con¬ 
sidered as rough estimates, we used very general and simple ap¬ 
proaches. We calculated the SFR based on the IR luminosity 
integrated between 8 pm and lOOOyum, 

SFR [Mo yr"‘] = 1.72 X 10"'° X 4 b , (2) 


following the relation for local galaxies from Kennicutt ( 1998| >. 

We estimated the black hole accretion rate (Mg^) from the 
equation 


Bol rlR _ yiVacc 2 
^AGN^AGN “ 


(3) 


based on the IR luminosity of the AGN. This relation as¬ 
sumes that a fraction of the rest-frame energy of the matter ac¬ 
creting onto the black hole is converted into radiation over the 
whole electromagnetic spectrum. The conversion is given by the 
efficiency factor 6, only slightly constrained by empirical stud¬ 
ies. While for example [Shankar et al.|(|2010|) suggest e > 0.2 
based on quasar clustering, others like[Davis & Faor[([2011[) find 
0.06 < 6 < 0.4, depending on the mass. Here, we used a con¬ 
servative number 6 = 0.1 since we were interested in upper lim¬ 
its for the black hole accretion rate. Furthermore, we needed to 
make an assumption for the factor /c, converting the IR AGN lu¬ 
minosity to the bolometric luminosity. This factor is uncertain 
and can vary between 1.4 and 15. We adopted the number from 
Drouart et al. ( |2014[ ), k = 6, assuming an RF AGN template. 

Figure [10] shows the estimated upper limit of the SFR and 
an estimate of the black hole accretion rate versus redshift based 
on the maximum IR luminosities derived in Sect. 15.2.21 For the 
accretion rate, we also showed the related uncertainty, arising 
from the assumptions made for the efficiency factor e and the 
conversion factor k as discussed above. Varying e between 0.06 
and 0.4, and k between 1.4 and 15, representing the range of 
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Fig. 10. Star formation rate (blue line) and accretion rate (red line) ver¬ 
sus redshift. The calcul ation of both rates is based on the modelling 
described in Sect. [5.2.2l matching the measured 1.4 GHz flux density of 
15 mJy and providing the lowest deviation to all available flux density 
upper limits. The red shaded area represents the uncertainty in the de¬ 
termination of the BH accretion rate, based on varied numbers for the 
efficiency factor e and the conversion factor k. We note that the SFR is 
meant to be an upper limit since the SB contribution is only constrained 
by upper limits. 


reasonable numbers, results in the uncertainties shown by the 
red shaded area in Fig.fTO 


We found that IFRS at z < 4.5 are constrained by SFRs 
of a few hundred solar masses per year. The SFR upper limit 
increases with redsh ift and exceeds lO'^Moyr"^ at z ~ 11. 
Drouart et al.[ ( |2Q14[ ) find SFRs for their HzRGs in the range of 


lOOMoyr"^ to 5000Moyr“k These numbers are in agreement 
with our results, considering that the SFRs calculated here are 
meant to be upper limits. 

The estimates of the black hole accretion rate are below 
1 Mq yr“^ at z < 3 and below 10 Mq yr“^ at z ~ 6. At redshifts 
z > 10 , the accretion rate is below a few tens of solar masses per 
year. Drouart et al. (2014) find accretion rates between 1 Mq yr“^ 
and lOOMoyr"^ for their sample of HzRGs. Again, our results 
are in agreement with the numbers presented by [Drouart et al. 


Considering that the HzRGs in the sample of [Drouart et al. 


are 

at z < 5, IFRS seem to show lower accretion rates if they are lo¬ 
cated at similar redshifts. In contrast to the SFR estimated in this 
work, which are meant to be upper limits as discussed above, the 
black hole accretion rates are considered to be rough estimates 
under the assumption that the SEDs of IFRS are similar to those 
of HzRGs. However, the assumption for the conversion factor k 
and the efficiency factor e and the scattering of the IR luminosity 
of AGNs ( [Drouart et al. |20 1 4} add significant uncertainties. 


All our findings in this IR SED modelling are in agreement 
with the results by [Drouart et al. ( 2^Q14[ ) for HzRGs. The redshift- 
dependent IR luminosity upper limits and the directly related 
SFRs and black hole accretion rates derived in this work agree 
with the numbers calculated by [Drouart et al.[ for their sources 
which, in contrast to IFRS, partially provide FIR detections. In 
our IR SED modelling, we did not find any evidence disproving 
the hypothesis that IFRS are similar to HzRGs. 
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6. Conclusion 

In this work, we presented the first FIR data of IFRS. Six IFRS 
have been observed between lOOyum and 500 yum with the in¬ 
struments PACS and SPIRE on board the Herschel Space Ob¬ 
servatory. None of the observed IFRS has been detected at any 
of the five FIR wavelengths down to median 3(T levels between 
4.4 mJy at 100yum and 17.6 mJy at 500yum. Even the stacking of 
the five maps of five IFRS at each wavelength did not provide a 
detection. 

We used the FIR flux density upper limits combined with the 
radio detections and the SERVS/SWIRE fiux density upper lim¬ 
its to model the broad-band SED of IFRS. The characteristics of 
IFRS can only be explained by known SEDs of HzRGs (e.g. the 
spiderweb galaxy) if these templates are shifted to z > 10.5 or 
by the SED of Cygnus A at 5 < z < 8.5. All other templates, for 
example RL HyLIRGs, ULIRGs with an AGN in their centre, or 
RL quasars fail to reproduce the SED constraints of IFRS. We 
also tested whether modified templates of known galaxies can 
explain the characteristics of IFRS at lower redshifts z < 5. In 
this approach, we scaled the templates in luminosity to match 
the measured fiux densities at 1.4 GHz or added extinction in the 
rest-frame optical and near-IR regime if required. We found that 
the templates of the spiderweb galaxy, Cygnus A, the CSS source 
3C48, and the RL quasar 3C273 can be modified to match all 
photometric constraints of IFRS. However, at low redshifts, ad¬ 
ditional obscuration by dust was needed for most SED templates 
to match the near-IR faintness of IFRS. 

Although no IFRS has been found at low redshifts z ~ 1, our 
modelling provided appropriate modified templates for this sce¬ 
nario. This implies that our modelling is physically unrealistic at 
low redshifts, potentially because of scaling down the luminosity 
and—by this—the SMBH mass which might prevent the object 
to be an RL AGN. On the other hand, if the IFRS analysed in 
this work are indeed at z ~ 1, they would form a new subclass 
of IFRS with different characteristics than the IFRS with known 
spectroscopic redshifts. Particularly, we showed that SED tem¬ 
plates of star forming galaxies, Seyfert galaxies, ULIRGs, and 
HyLIRGs are inappropriate to reproduce the photometric con¬ 
straints of IFRS. 

We modelled the maximum IR SED of IFRS based on a set 
of SB templates and an AGN contribution and using all available 
fiux density upper limits, aiming at measuring the maximum IR 
fiux of IFRS. We found that this model is unphysical at lower 
redshifts because of its disagreement with the radio-IR corre¬ 
lation. We found that the IFRS could also be modelled by the 
spiderweb galaxy template, together with an additional SB com¬ 
ponent. The related IR luminosity upper limits are in agreement 
with those of HzRGs. Using these maximum IR luminosities, 
we estimated black hole accretion rates and upper limits for the 
SFR. These numbers agree with the numbers of HzRGs, too. 

In summary, the IFRS discussed here might be (a) objects 
identical to known HzRGs, but at high redshifts (z > 10.5), 

(b) objects similar to Cygnus A, but at high redshifts 5 < z < 8.5, 

(c) objects similar to known HzRGs or Cygnus A, but scaled 
down in luminosity, or (d) objects similar to CSS sources or RL 
quasars, but modified by significant additional extinction and 
scaled in luminosity which makes these templates to HzRGs, 
too. In any case, objects which reproduce the characteristics 
of IFRS fulfil the selection criteria of HzRGs. We estimated 
that IFRS contain SMBHs accreting at the rate of 1 Moyr“^ to 
50Moyr“^ together with star formation at a rate of up to sev¬ 
eral thousand solar masses per year in case of an additional SB 
component. 
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